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ABSTRACT 



Recent studies have shown that the elemental abundances in the Sun are anomalous when compared to most (about 85 %) nearby 
solar twin stars. Compared to its twins, the Sun exhibits a deficiency of refractory elements (those with condensation temperatures 
7*c <; 900 K) relative to volatiles (T c ^ 900 K). This finding is speculated to be a signature of the planet formation that occurred more 
efficiently around the Sun compared with the majority of solar twins. Furthermore, within this scenario, it seems more likely that the 
abundance patterns found are specifically related to the formation of terrestrial planets. In this work we analyze abundance results from 
six large independent stellar abundance surveys to determine whether they confirm or reject this observational finding. We show that 
the elemental abundances derived for solar analogs in these six studies are consistent with the T c trend suggested as a planet formation 
signature. The same conclusion is reached when those results are averaged heterogeneously. We also investigate the dependency of 
the abundances with first ionization potential (FIP), which correlates well with T c . A trend with FIP would suggest a different origin 
for the abundance patterns found, but we show that the correlation with T c is statistically more significant. We encourage similar 
investigations of metal-rich solar analogs and late F-type dwarf stars, for which the hypothesis of a planet formation signature in the 
elemental abundances makes very specific predictions. Finally, we examine a recent paper that claims that the abundance patterns 
of two stars hosting super-Earth like planets contradict the planet formation signature hypothesis. Instead, we find that the chemical 
compositions of these two stars are fully compatible with our hypothesis. 

Key words, stars: abundances - Sun: abundances - stars: planetary systems 



1. Introduction 

Finding planets outside the solar system is one of the major en- 
deavors of contemporary astrophysics. Large efforts to detect 
them using radial velocity, transits, and microlensing observa- 
tions are currently underway as they have been proven success- 
ful. Indeed, more than 450 extrasolar planets have been discov- 
ered to date. A majority of the planets that have been detected 
have masses similar to that of Jupiter and a significant fraction 
of them are in close-in, short-period orbits, hen ce their classifi- 
cation as "hot Jupiters" (see, e.g., the review by|Ud rv & S antos 
|2007|) . The parent stars are mostly G and K dwarfs. Due to the 
intrinsic limitations of the techniques currently used to find exo- 
planets, however, terrestrial planets in Earth-like orbits have re- 
mained elusive. This situation may change in the next few years 
as transiting data from the Kepler mission (iBasri et al.ll2005l) are 
scrutinized spectroscopically. In the meantime, we face the chal- 
lenge of developing new practical techniques to find those ob- 
jects. 

Comparisons of the chemical composition of stars that host 
extrasolar planets and stars for which planets have not yet 
been detected have been performed by a number of groups 
(e.g., Gonzalez et al. 2001; Sadakane et al. 2002; Heiter & Luckl 
2003; iSantos et al.l2004UFischer & Valentil2005l) . These studies 
unanimously agree on the metal-rich nature of planet hosts (i.e., 
stars known to host planets are on average more metal-rich than 



a typical field star), first suggested bv lGonzalezI (1 19971) . From a 
practical point of view, this result can be used to increase the ef- 
ficiency of exoplanet searches. Note, however, that it is not clear 
whether evolved stars which host planets ar e, on average, more 
metal-rich than a typical field star (e.g., Hekker & Melendez] 
120071: iPasquini et alj|2007tlTakeda et alj|2008l) . 

Attempts have been made to find other indicators in the pho- 
tospheric chemical composition of stars that would suggest the 
presence of extrasolar planets. For instance, it has been sug- 
gested that the lithium abundances in stars tha t host planets 
are more depleted than those in single stars ( e.g., IChen & Zhao! 
2006; iGonzalezI 120081: llsraelian et ail 120091) as a result of a 
star-planet interaction that has led to enhanced mixing and/or 
deepening of the convective zone. This, in turn, increases the 
rate of lithium burning in the stellar interior. It is possible, 
however, that these results have been affected by a combi- 
nation of heterogeneous abundance results or, more impor- 
tantly, observational biases related to s tellar age and tnetal- 
licitv (iRvanl 120001; iLuck & Heiterll2006l IMelendez et all 120101: 



Bau mann et al.l 2*010l) . Therefore, the possible connection be- 



tween lithium abundance and presence of planets is not conclu- 
sive. 

IMelendez et al.l (l2009i hereafter lM09l) have shown that the 
chemical composition of the Sun is anomalous when compared 
to stars of very similar fundamental parameters, so called "solar 
twins." Inspection of abundance trends with condensation tem- 
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perature (7c) revealed, for the first time, that, compared to most 
solar twins (about 85 %), the Sun has a deficiency of refractory 
elements (T c £ 900 K) relative to volatiles (T c < 900 K). In 
the solar photosphere, elements with the highest Tc are depleted 
by about 20 % relative to volatiles. A 20 % difference in relative 
abundances (about 0.08 dex in the standard [X/H] scalefl is ex- 
tremely difficult to detect because systematic errors are often of 
similar magnitude or larger. In the M09 study, those errors were 
minimized by analyzing only stars with fundamental parameters 
very similar to solar, for which systematic errors in the derived 
abundances cancel-out in a strictly differential analysis. Using 
a different sample of stars and a new set of observations (solar 
twins in the northern hemisphere as opposed to south ern hemi- 
sphere stars as in the M0^case),|Ramfrez et al. I J20091 hereafter 
R09) determined abundance trends with condensation te mpera- 
ture which are in excellent agreement with those found by M09. 
Both the amplitude of the effect (^ 20 % deficiency of refracto- 
ries relative to volatiles) and the frequency of stars in which the 
trend is obser ved (^ 85 %) were con firmed by this latter s tudy- 
Mo re rec ently, iGonzalez et at] d2010l) have confirmed both lM09l 
and |R09| results using a different sample of stars and approach 
to the analysis of the stellar abundance data. In addition, they 
explore the impact of effective temperature and metallicity on 
these abundance trends. 



Table 1. Samples analyzed in this paper, selection criteria for 
solar analogs, and number of stars selected. 



lM09l have speculated that this "solar abundance anomaly" 
is due to the formation of planets in the solar system, and pos- 
sibly specifically the terrestrial planets, because the missing re- 
fractories are currently inside them. If this hypothesis is proven 
correct, very precise relative abundances could be used to iden- 
tify stars that harbor terrestrial planets. A lack of correlation in 
the abundances relative to solar with Tc would indicate the pres- 
ence of such systems around other stars. The amount of observ- 
ing time required to identify planet hosts with this "abundance 
method" is orders of magnitude smaller than that necessary to 
detect planets using radial velocity or transit techniques. Only 
one set of observations is required per star, in contrast to the 
lengthy campaigns that have to be carefully designed to detect 
very small radial velocity or photometric variations. 

The identification of terrestrial planet hosts using abun- 
dances will help us to answer important questions such as (1) 
how rare is the solar system? and (2) which types of stars are 
more likely to host terrestrial planets? However, the connection 
between the abundance trends with condensation temperature 
and terrestrial planet formation needs to be explored further. 

A large number of stellar photospheric chemical composi- 
tion studies have been published to date. Only a few, however, 
are based on a homogeneous analysis of high quality data for 
large samples of stars and deal with more than a handful of 
chemical elements. In this paper, we compile results from six 
of those studies. We concentrate on stars that are similar to 
the Sun, since the signature is difficult to detect and possibly 
blurred by systematic errors in larger samples including stars 
with significantly different fundamental parameters. We analyze 
the trends with condensation temperature for these samples and 
check whether other independent st udies of s tellar abundances 
support the observational findings of M09 and R09. 



Sample 


Solar analogs range 


Number 




A7W/Alogg/A[Fe/H] 


of stars 


Reddv et al. (2003) 


+300 K/ +0.2/ +0.2 


11 


Allende Prieto et al. (2004) 


±150K/±0.2/+0.2 


12 


Takeda (2007) 


±150K/±0.2/+0.2 


14 


Neves et al. (2009) 


+ 100K/+0.1/+0.1 


19 


Gonzalez et al. (2010) 


+ 100K/+0.1/+0.1 


11 


Bensbv et al. (in preparation) 


±100K/±0.1/+0.1 


20 



2. Sample selection and abundances adopted 

We ha ve comp iled abundance res ults f r om six independent 
studie s: iReddv et all (120031 hereafte r iROll). I Allende Prieto et al.l 
(I2004L her e after lA04h iTakedal 120071 hereafte r |T07l) 
lNevesetal.1 (120091 hereafter lN09h IGonzalez et al.l (2010, 
herea fter iGlOl) . and iBensbv et alj (in preparation; hereafter 
BFO). A common characteristic among these studies is the 
homogeneity of the analysis and the high quality of the data 
employed. Each has adopted a unique temperature scale, 
spectral line selection, and method of abundance determination. 
In addition, abundances were determined from equivalent width 
matching and/or line-profile fitting using spectra of high resolu- 
tion (R > 50,000) and high signal-to-noise ratio (R > 100). In 
all cases, abundances for at least 10 elements were determined. 

From each sample, only those stars with fundamental param- 
eters (r e ff, logg, [Fe/H]) within a certain range from the solar 
values (Tig = 5777K,logg = 4.44, [Fe/H] = 0) were selected. 
In the R09 study, the ranges adopted were ±100K in effective 
temperature, +0.1 dex in surface gravity, and ±0.1 dex in iron 
abundance. An identical selection was made for the N09,|GlS 
and|BFO samples. For the R03~l lA04l and lT07l samples, however, 
we used wider ranges of allowed stellar parameters in order to 
have at least 10 solar analogs, as detailed in Table [TJ The num- 
ber of stars selected from each study are listed also in Table[TJ In 
Fig. Q] we illustrate the overall properties of each sample, high- 
lighting the location of the stars selected for this work. 



We use the standard notation: Ax = lognx/«H+12, where nx is the 
number density of the element X ; [X/H] = A x - A°; and [X/Fe] = 
[X/H] - [Fe/H]. 



R03 used spectra obtained with the Robert G. Tull spectro- 
graph on the 2.7 m Telescope at McDonald Observatory. The 
main goal of this work was to study Galactic chemical evolu- 
tion using 181 nearby stars, almost all of them members of the 
Galactic thin disk. They derived abundances for 27 elements. 
The|R03 results were complemented with data for 1 76 stars with 
high p robability of being thick-disk members by Redd v~et~al1 
(2006). If we include the data from this later study in our work, 
the precision of the average abundance ratios that we derive for 
solar analogs (cf. Sect. [3) decreases from about 0.05 to 0.09 dex. 
The abundance trend we find with condensation temperature be- 
comes more noisy but i ts shape does not change significantly. 
We did not include the IReddv et a l. (2006) results in our work 
because we require the highest precision possible. Small changes 
in the determinati on of stellar parameters between IroI and 
IReddv et alJ (H006) and the slightly lower quality of the data em- 
ployed by the latter may be responsible for the increased scatter 
of the combined sample. 

A04 performed a volume limited survey of 1 1 8 stars more lu- 
minous than My = 6.5 mag within a 15 pc radius from the Sun. 
Their data were obtained at McDonald Observatory and wit h the 
FEROS spectrograph on the 1 .52 m Telescope at La Silla. 
derived abundances for 16 elements. The spectra from the 
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Fig. 1. Upper panels: effective temperature (Teg) versus surface gravity (logg) for all stars in the IroI IXoH IN09L iGlOl and|BFQ| 
studies (crosses). Filled circles represent the solar analogs used in this work. Lower panels: T e g versus [Fe/H] for the stars shown in 
the upper panels. 
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survey are available onlineQ We re-determined the stellar pa- 
rameters r e »f,logg, [Fe/H] of a number of solar analogs using 
these data and the methods described in R09. We found a mean 
difference of 152 K in T e g and 0.06 dex in log g between the re- 
derived parameters and those given by|A0S with the latter being 
systematically s maller . Indeed, the zero point of the temperature 
scale adopted bv lA04l has been recently revise d and found to be 
signif i cantly cooler than previously thought (Casagrande et al. 
120101) . |A04| provide detailed tables that allowed us to compute 
corrections to the abundances of all elements given an effective 
temperature correction. Thus, we applied a 4152 K correction to 
the parameters and abundances from IA041 We did not correct 
the abundances for the shift in log g values because there are no 
corresponding correction tables. The selection of solar analogs 
from this study was made after the r e ff correction to the abun- 
dances was applied. This type of correction was not possible for 
the other samples because the observational data are not publicly 
available, thus preventing us from determining the temperature 
correction required, and because none of the other studies pro- 
vide detailed tables of abundance corrections for given shifts in 
fundamental parameters. We discuss the implication of this cor- 
rection for the lA04l results to our work in Sect. [3] 

Fundamental parameters and elemental abundances of 15 el- 
ements for 160 FGK dwarf and subgiant stars observed with the 
1.88 m Telescope at Okayama Astrophysical Observatory were 
presented bv lT07l This work discusses chemical evolution, com- 
parison of abundances in planet-host and single stars, and empir- 
ical tests of the validity of the LTE approximation in abundance 
studies. 

Using spectra from the HARPS/GTO planet search pro- 
grams, [N09| determined abundances of 12 elements in 451 stars. 
The stellar parameters T s g, log g, [Fe/H] adopted in this study 



2 Spectroscopic Survey of Stars in the Solar Neighborhood (S 4 N): 
http://leda.as.utexas.edu/s4n/ 



are from.Sou sa et al . (2008). The wavelength coverage and typ- 
ical S/N of the HARPS data is unfortunately not appropriate to 
derive precise abundances of the volatile elements C, N, and O. 
However, ve ry pre cise abundances for many refractory elements 
are given bvlN09l 

For more than a decade, G. Gonzalez and collaborators 
have collected high quality spectra and analyzed them homo- 
geneously to derive precise stellar parameters and abundances 
with the goal of investigating the con nection between stellar 
abundances and extr asolar planets (e.g.. lGonzalezlll99"7l 1998; 
iGonzalez et al.ll200lT) . In their most recent publication ( GlOh . 
the abundance correlation with condensation temperature is ex- 
plored for about 85 % of their 85 stars with planets and 59 stars 
not known to host planets. The approach used in IGlOl consists 
of a differential comparison of abundance trends weighted by 
th e distance betwee n stars in stellar parameter space, as defined 
in iGonzaiezl (T2008 ) . In this way they are able to perform differ- 
ential analyses not restricted to solar twins. Abundances for ten 
elements are available from lGlOl although a comment must be 
made about their Sc abundances, which are based on two Sc 11 
lines. In lR09l three Sci lines were p referr ed. The reason was 
that the star-to-star scatter obtained by R09 for their solar twins 
reduced from 0.038 dex for Sc n to 0.018 dex for Sc 1. The aver- 
age [Sc/Fe] abundance ratio for Sci lines was 0.024 dex while 
that for Sc 11 lines was -0.003. Although non-LTE effects are ex- 
pected to be stronger for the Sc 1 lines, the [Sc i/Fe] ratio, which 
is the quantity that we use in this paper, might be less sensitive 
to non-LTE because the iron abundances are dominated by the 
large number of Fe 1 lines. It is likely that the systematic errors 
are more similar between Fe 1 and Sc 1 lines than between Fe 1 and 
Sc 11, making Sc 1 lines more reliable for strictly differential anal- 
ysis. We proceed our study including the Sc abundances from 
iGlOl but we will also comment on how their exclusion, justified 
by the statements made above, affects the conclusions. 
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iBensbv et al.1 (120031 120051) presented abundances for 12 el- 
ements in 102 F and G dwarf stars in the solar neighborhood 
based on high-resolution spectra obtained with the FEROS/ESO 
1.5m, SOFIN/NOT, and UVES/VLT spectrographs. Expanding 
those studies, the lBFOi study will contain a total of 703 F and G 
dwarf stars, adding 600 new stars that have been observed with 
the MIKE spectrograph on the 6.5 m Clay Magellan telescope. 
These new spectra cover the full optical wavelength regime 
(350-920 nm), have high spectral resolution (R a 65,000), and 
hi gh sig nal-to-noise ratio ( S IN > 250). As the a nalys is method 
in BFO is very similar to IBensbv et al.1 d2003L 120051) . the un- 
certainties in stellar parameters and abundance ratios are simi- 
lar as well. First results from this e xpanded data set have bee n 
pre sented inlBensbv et alj d2007allbh . iFeltzin g & Bensbvl d2008l) . 
and B ensbv & Feltzingl(l2010l). The relev ant atomic data adopted 
bv lBFOl are given in lBensbv et al . (2003). 

The oxygen abundances from the iBFOl work were derived 
by them from the 777 nm O i triplet line s, which are know n 
to be severely affected by non-LTE (e.g.. Fab bian et al.l l2009). 
Although the range in stellar parameters adopted for the solar 
analogs in this sample is small, |R09| have found that differen- 
tial non-LTE corrections slightly improve the results. Therefore, 
we corrected the oxygen abundances froml BFO using the non- 
LTE correction tables bv lRamfrez et al.l ( 12007b . We find that the 
star-to-star scatter in the average [O/Fe] abundance ratio that we 
derive for the solar analogs from the lBFOi study decreases from 
0.048 dex in LTE to 0.038 dex in non -LTE. We did not use the 
more recent non-LTE corrections by Fabb ian et alJ (2009) be- 
cause their tables do not cover [Fe/H] > 0. The differential n on- 
LTE correction (star relative to Sun) bv lRamfrez et al. I (12001 for 
a star that is 100K warmer than the Sun (the ot her parameters 
being solar) is about 0.01 dex greater than that by Fabbian et al. 
(2009). A similar calculation for a [Fe/H] = -0.1 star yields a 
difference of -0.02 dex. The impact of logg differences is neg- 
ligible within the range allowed for solar analogs in this paper. 
These differences are small but could have an impact on the av- 
erage [O/Fe] abundance ratios we derive. Unfortunately, we are 
unable to use lFabbian et al.l ([2009 ) non-LTE correction tables for 
the reason expos ed above. 

As shown by Bensbv et al. ( 2003, 2005), abundance ratios of 
kinematically selected thin- and thick-disk stars are very similar 
at [Fe/H] ^ 0.0. This is not the case at lower metallicities. We 
checked that the results presented in Sect. [3] are not affected by 
the inclusion of thick-disk solar analogs (i.e., excluding them 
did not have a significant impact). Only about 15 % of the stars 
studied in this paper have high kinematic probability of being 
thick-disk members, a number that is roughly consistent with 
current estimates of the fraction of thick-disk stars in the solar 
neighborhood. 

3. Abundances versus condensation temperature 

For each of the six samples analyzed in this paper we deter- 
mined an average [X/Fe] vs. Tq relation. For each element X, 
we calculated an average [X/Fe] and standard deviation with the 
data for the solar analogs (weighted by the observational errors, 
when available). We adopted the 50 % equilibrium condensation 
temperatures fo r a solar system composition gas computed by 
lLoddersI (120031) . The resulting average [X/Fe] vs. Tq relations 
ar e show n in Fig. [2] along with the corresponding relation found 
bv lR09l Table|2]lists the average [X/Fe] abundance ratios of so- 
lar analogs that we obtain for each of the six studies used in this 
paper. We also list in Table |2]our adopted condensation temper- 
ature and first ionization potential for each element. The latter 



are adopted from the "Atomic Properties of the Elements" com- 
pilation by the National Institute of Standards and Technology 
(NIST)0 In Table [3] we show the average [X/Fe] abundance 
ratios of solar twin stars from the MQ9| and lR09l studies. The 
adopted condensation temperatures and first ionization poten- 
tials for the elements studied in these two papers are given in 
Table 12 except for element P (Z=15), which is only covered by 
the lM09l study (T c = 1229 K, FIP = 10.49 eV). 

The abundances of volatile elements in the solar analogs 
from the R03 study are in severe disagreement with those found 
bv lR09l but their Tq trends for the refractory elements are in rea- 
sonably good agreement (Fig.|2^). In this study, the spectroscop- 
ically determined abundances, relative to iron, of carbon, nitro- 
gen, and oxygen, which are the elements with the lowest Tq in 
Fig.[2k, are very sensitive to errors in the T e g scale adopted. The 
T e g sensitivity of these elements is stronger than that of other 
elements. R03 noticed, for example, that their [C/Fe] abundance 
rat ios are systematically hi gher by 0.15 dex compared to those 
by lEdvardsson et al.l ( 1 1993b . and attributed the difference to an 
offset of about 80 K in the T e g scales adopted. Their [N/Fe] abun- 
dance ratios were determined from the analysis of two weak N 
lines. Errors in the continuum determination of weak lines may 
lead to severely overestimated equivalent widths when noisy 
data are used. In these cases, relative abundances can be system- 
atically off if the reference spectrum (the solar one) has much 
higher signal-to-noise ratio, as is often the case. Indeed, lR03l 
noted that their N abundances are abo ut 0.2 dex h i gher o n av- 
erage compared to those determined bv lClegg et al.l ( 198lT) . who 
used also atomic N lines. The other two elements that depart 
significantly from the lR09l result are Zn and S. In addition to po- 
tential systematic and statistical errors which could have re sulted 
in R03 abundances being significantly higher than those bv lR09l 
it should be noted that this sample is, on average, slightly more 
metal-poor than the Sun (Fig. [TJ. This behavior is not observed 
in any of the other samples of solar analogs, which are shown 
to distribute uniformly around [Fe/H] = 0. Thus, it is also pos- 
sible that small chemical evolution effects are partially respon- 
sible for the discr epanc y observed. Indeed, the four lowest Tc 
elements from the |R03| work show decreasing abundance ratios 
with [Fe/H], although the s lopes are dissimilar for the different 
elements (see Fig. 9 in lR03l) . 

The A04 abundances seem to broadly agree with the solar 
abundance anomaly reported by M09 and IR09L despite the rel- 
atively large uncertainties (Fig. %2&). On average, volatiles have 
[X/Fe] ^ -0.06 whereas refractories have an average [X/Fe] 
of about 0.02 (weighted average). The Tc trend appears to be 
steeper compared to R09 but the overall behavior of the abun- 
dances with Tc is qualitatively similar. Since we shifted the pa- 
rameters and abundances from |A04| following a 152 K offset 
in T e g (Sect. ©, it is important to investigate the results using 
the original values. In this case the average abundance ratio of 
volatiles is about 0.08 dex whereas that corresponding to the re- 
fractories is about 0.09 dex. The linear Tc trend for the refrac- 
tories is conserved but the difference in average abundance be- 
tween volatiles and refractories decreases significantly. Th e re- 
sult using the original abundances resembles somewhat the lR03l 
result (although the abundances of volatiles are not that high 
in the |A04| case), which reinforces our suspicion that the abun- 
dances of volatiles in these two studies (the original ones in the 
lA04l case) are too high due to errors in the T e g scale. As men- 
tioned in Sect. [2l the T e g scale a dopted in these two studies, 
namely that by lAlonso et alJ (119961) . has been recently shown 



3 Available online at http://www.nist.gov/physlab/data/periodic.cfm 
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Table 2. Average (star-to-star) elemental abundance ratios [X/Fe] of solar analogs in each of the samples studied in this work. Error 
bars are the weighted 1-cr scatter. 



z 


El. 


T c 


FIP 


R03 


A04 


T07 


N09 


G10 


BFO 






(K) 


(eV) 






6 


C 


40 


11.26 


0.138 ± 


0.034 


-0.044 + 0.111 










7 


N 


123 


14.53 


0.238 ± 


0.052 












8 


O 


180 


13.62 


0.048 ± 


0.046 


-0.132 ±0.075 








-0.008 ± 0.038 


11 


Na 


958 


5.14 


-0.011 ± 


0.051 




-0.053 ± 0.054 


-0.032 ± 0.037 


-0.067 ± 0.030 


-0.027 ± 0.041 


12 


Mg 


1336 


7.65 


0.005 ± 


0.052 


0.002 ± 0.063 


-0.012 ±0.051 


0.035 ± 0.034 




-0.028 ± 0.035 


13 


Al 


1653 


5.99 


0.035 ± 


0.041 




-0.008 ± 0.055 


0.031 ± 0.040 


-0.026 ± 0.046 


-0.020 ± 0.042 


14 


Si 


1310 


8.15 


0.031 ± 


0.030 


0.010 ± 0.030 


0.005 ± 0.041 


0.011 ±0.030 


-0.039 ± 0.030 


-0.018 ±0.030 


16 


S 


664 


10.36 


0.148 ± 


0.037 




-0.053 ± 0.061 








19 


K 


1006 


4.34 


0.070 ± 


0.074 












20 


Ca 


1517 


6.11 


-0.016 ± 


0.030 


-0.014 ± 0.070 


0.006 ± 0.030 


0.033 ± 0.030 


-0.001 ± 0.030 


0.016 ± 0.030 


21 


Sc 


1659 


6.54 


0.047 ± 


0.066 


0.073 ± 0.048 


0.042 ± 0.080 


0.028 ± 0.030 


-0.061 ± 0.041 




22 


Ti 


1582 


6.82 


-0.025 ± 


0.030 


0.097 ± 0.047 


0.017 ±0.030 


0.039 ± 0.030 


-0.007 ± 0.030 


0.004 ± 0.030 


23 


V 


1429 


6.74 


-0.067 ± 


0.030 




-0.014 ± 0.054 


0.014 ± 0.030 


-0.020 ± 0.030 




24 


Cr 


1296 


6.77 


-0.020 ± 


0.030 




0.013 ± 0.030 


0.006 ± 0.030 


-0.005 ± 0.030 


-0.013 ± 0.030 


25 


Mn 


1158 


7.43 


-0.085 ± 


0.039 




-0.015 ± 0.056 


-0.023 ±0.031 




. . . 


27 


Co 


1352 


7.86 


-0.047 ± 


0.051 


-0.012 ± 0.033 


-0.010 ± 0.050 


0.001 ± 0.030 


-0.041 ± 0.042 




28 


Ni 


1353 


7.63 


-0.023 ± 


0.030 


-0.007 ± 0.030 


-0.008 ±0.031 


-0.010 ±0.030 


-0.035 ± 0.030 


-0.020 ±0.031 


29 


Cu 


1037 


7.73 


-0.070 ± 


0.067 


-0.032 ± 0.062 


-0.058 ± 0.063 








30 


Zn 


726 


9.39 


0.055 ± 


0.042 


0.040 ± 0.089 


-0.005 ± 0.055 






-0.018 ± 0.052 


38 


Sr 


1464 


5.70 


-0.008 ± 


0.080 












39 


Y 


1659 


6.38 


0.077 ± 


0.060 


0.113 ±0.090 








0.010 ± 0.047 


40 


Zr 


1741 


6.84 


0.047 ± 


0.118 












56 


Ba 


1455 


5.21 


0.050 ± 


0.032 


0.095 ±0.150 








0.066 ± 0.080 


58 


Ce 


1478 


5.47 


0.055 ± 


0.071 


0.156 ±0.100 










60 


Nd 


1602 


5.49 


0.085 ± 


0.057 


0.133 ±0.090 










63 


Eu 


1356 


5.67 


0.040 ± 


0.051 


0.086 ± 0.056 











Table 3. Average (star-to-star) elemental ab undance ratios 
rX/Fel (and 1-c r scatter) of sol ar twins in the iRamfrez et al.l 
(2009,|R0l) and lMelendezltaD (12009. M09) studies. 



z 


El. 


M09 


R09 


6 


C 


-0.048 


± 0.042 


-0.065 ± 


0.033 


7 


N 


-0.060 


±0.051 






8 


O 


-0.033 


± 0.038 


-0.015 ± 


0.028 


11 


Na 


-0.024 


± 0.024 


-0.052 ± 


0.027 


12 


Mg 


-0.018 


± 0.030 






13 


Al 


0.033 


± 0.025 


-0.012 ± 


0.053 


14 


Si 


-0.030 


± 0.020 


-0.012 ± 


0.019 


15 


P 


-0.051 


± 0.063 






16 


S 


-0.056 


± 0.042 


-0.036 ± 


0.028 


19 


K 


-0.013 


± 0.020 






20 


Ca 


-0.004 


±0.010 


0.024 ± 


0.019 


21 


Sc 


0.018 


±0.031 


0.024 ± 


0.018 


22 


Ti 


0.011 


±0.016 


0.028 ± 


0.031 


23 


V 


-0.009 


±0.019 


0.014 ± 


0.015 


24 


Cr 


-0.011 


±0.012 


0.016 ± 


0.017 


25 


Mn 


-0.029 


±0.018 


-0.018 ± 


0.021 


27 


Co 


-0.011 


± 0.028 






28 


Ni 


-0.007 


±0.013 


-0.021 ± 


0.013 


29 


Cu 


-0.021 


± 0.023 


-0.025 ± 


0.033 


30 


Zn 


-0.023 


± 0.032 


-0.027 ± 


0.061 


39 


Y 


0.015 


± 0.033 


0.024 ± 


0.043 


40 


Zr 


0.023 


± 0.027 


0.037 ± 


0.060 


56 


Ba 


0.012 


± 0.025 


0.070 ± 


0.055 



to suffer fr om a relatively large sys tematic error (^ 100 K) in its 
zero point (Casagran de et al.ll2010h . 

Fig. [2J: shows that the abundance data for solar analogs by 
|T07| are consistent with the [X/Fe] vs. Tc trend found in our 
previous papers. The most volatile elements (C, N, O) are not 



Table 4. Average (star-to-star) elemental abundance ratios (and 
1-cr scatter) of solar analogs from the six independent studies 
examined in this work. The second column gives the average 
abundance ratios for the purely spectroscopic papers only. 



z 


El. 


average (all) 


average (spec.) 


6 


C 


-0.044 ±0.111 






8 


O 


-0.033 ±0.071 


-0.008 


± 0.038 


11 


Na 


-0.048 ± 0.020 


-0.048 


± 0.020 


12 


Mg 


0.002 ± 0.026 


0.001 


± 0.029 


13 


Al 


0.005 ± 0.027 


-0.003 


± 0.024 


14 


Si 


-0.000 ± 0.023 


-0.012 


± 0.020 


16 


S 


-0.053 ± 0.061 


-0.053 


± 0.061 


20 


Ca 


0.007 ± 0.020 


0.014 


± 0.020 


21 


Sc 


0.042 ± 0.020 


0.030 


± 0.020 


22 


Ti 


0.012 ±0.032 


0.013 


± 0.020 


23 


V 


-0.023 ± 0.032 


-0.004 


± 0.020 


24 


Cr 


-0.004 ± 0.020 


0.000 


± 0.020 


25 


Mn 


-0.042 ± 0.030 


-0.021 


± 0.020 


27 


Co 


-0.016 ±0.020 


-0.013 


± 0.020 


28 


Ni 


-0.017 ± 0.020 


-0.018 


± 0.020 


29 


Cu 


-0.045 ± 0.020 


-0.058 


± 0.063 


30 


Zn 


-0.004 ± 0.020 


-0.012 


± 0.020 


38 


Sr 


-0.008 ± 0.080 






39 


Y 


0.047 ± 0.041 


0.010 


± 0.047 


40 


Zr 


0.047 ±0.118 






56 


Ba 


0.054 ± 0.020 


0.066 


± 0.080 


58 


Ce 


0.089 ± 0.067 






60 


Nd 


0.099 ±0.031 






63 


Eu 


0.061 ±0.032 







covered by lT07l but the break in the Tq trend at about Tc — 
1000K, which is crucial for the interpretation of this abundance 
trend (Sect.|4]i, appears to be detected. 



6 



Ramirez et al.: A signature of terrestrial planet formation? 




500 



1000 
T (K) 



1500 



Fig. 2. Abundance ratios relative to iron for the solar twins and 
analogs from the six studies used in this paper as a function of 
condensation temperature (open circles). Filled squares with er- 
ror bars represent the average value and standard deviation for 
each element. Solid lines represent the average [X/Fe] value of 
volatiles (7c < 900 K) and the linear fit to the [X/Fe] versus 7c 
relation of refractori es (Tc > 900 K) derived for solar twins by 
iRarmrez et all d2009t) . 



The abundance vs. Tc trend obtained for the N09 sample 
of solar twins is in excellent agreement with IR09I and M09 
(Fig. |2jl). Note, however, that this comparison is limit ed to re- 
fractory elements since no volatiles were analyzed by NOjJ In 
any case, both the overall effect detected by |M09| and lR09l and 
the value of the [X/Fe] vs. T c slope for Tc > 900 K are in 
good agreement with the N09 data for solar analogs. Moreover, 

2 - 4 of the 19 solar analogs from |N09| show a [X/Fe] 
vs. Tc slope close to zero (the uncertainties in the abundance 
ratios of individual stars prevent us from determining this num- 
ber more accurately). This frequency of stars which do not fol- 
low the Tc trend (10-20%) is also in good agreement with the 
^ 15 % frequency found by both IR09I and M09. A previous in- 
vestigatio n of Tc trends by th e same group has been published 
before dEcuviUon et al] 12006). This previous work was unable 
to reach strong conclusions regarding Tc trends, probably be- 
cause no investigation of solar twin stars alone was made. This 
study also adopted s olar abundances from other works, mainly 
Ande rs & Grevessel (Il989l) . instead of performing a strictly dif- 
ferential analysis using their own solar spectra. We have shown 
that both the sample restriction to solar twins and strictly dif- 
ferential analysis are necessary to detect the small Tc trend in 
the abundances. T he fact that this group (Ecuv illon et al] [2006; 
iNeves et al. 2009) has now been able to detect the abundance 
signature using our approach ( Gonzale z Hernandez et al 1 120101) 
confirms this statement. 

In the Tc range between about 900 and 1600K, the agree- 
ment between the lR09l [X/Fe] vs. Tc trend and that given by the 
iGlOl solar twin data is good, albeit a very small offset in [Fe/H] 
of at most 0.02 dex (Fig. [2^). The two elements with the highest 
7c i n the l GlOl work are Al and Sc. The aver age A l abundance 
from IGlOl is in good agreement with that by lR09l and they are 
both within about 1-cr of the mean [X/ Fe] v s. Tc trend. This 
is not the case for the Sc abundance by IGIOL As explained in 
Sect. [2] this discrepancy may be due to the use of Sc n lines in 
the IGlOl work as opposed to Sci lines as in lR09l The latter were 
shown to reveal less star-to-star scatter for the solar twins and 
therefore they are expected to be more reliable in strictly differ- 
ential analysis. With the exception of the Sc a bunda nce s, the Tc 
trend in the IGlOl data is consistent with that in lR09i and lM09L 

Although the average [X/Fe] v s. Tc trend obtained for the 

BFO data is in agreement with the lR09l result within 1-cr, it is 
cle ar that the Tc trend of refractories is slightly more shallow in 
the lBFOl case (Fig. |2f). The difference in the mean abundance of 
refractory and volatile elements, weighted by the errors, is only 
about 0.01 dex, but the sign of such difference is in agreement 
with all other studies presented thus far in this paper as well as 
lM09l andlR09l 

In Fig. [3] w e show an average [X / Fe] vs. 7c relation for 
all six samples dROl lA04t lT07t lN09b IGlOl iBFOl) combined. 
For the elements present in more than one study we computed 
the weighted mean of the individually averaged results. We ex- 
cluded the lR03l results for the volatiles because they are likely 
affected by systematic errors in the temperature scale, measure- 
ment of equivalent widths, and/or chemical evolution effects, as 
explained above. We also excluded the K abundances from this 
study because of their large internal error (=; 0.07 dex) and the 
fact that they are based on the analysis of only one very strong 
feature (EW 170 mA), which could be affected by continuum 
placement and saturation. The impact of the mean Sc abundance 
from the IGlOl solar analogs, which we suspect to be slightly off, 
on this result is minor but we excluded it for consist ency. The 
agre ement of this average [X/Fe] vs. Tc trend with the lM09l and 
R09 results is good, but we emphasize the fact that this is an 
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Fig. 3. Abundance ratios relative to iron as a function of con- 
densation temperature for all samples of solar twins and analogs 
investigated in this paper (weighted average and l-cr scatter of 
mean abundance ratios from different studies). Solid lines rep- 
resent the a verage TX/Fel versus Tc relation for solar twin stars 
obtained by Ramire zetal] d2009l) . 




-0.15 



Fig. 4. As in Fig. [3] for the four more recent papers, which use 
only the spectra (and not photometry and/or trigonometric par- 
allaxes) to determine stellar parameters. 



heterogeneous compilation of abundances and that the different 
samples used for this compilation exhibit a slightly different be- 
havior when analyzed separately. In fact, we argue that because 
the data used to make this figure are heterogeneous and still show 
a clear a bu ndanc e vs. 7c correlation, it strengthens the findings 
of M09 and R09. The [X/Fe] abundance ratios for the average of 
the six samples combined are given in Table [4] 

Ins pectio n o f Fig . El clearly reveals that the results obtained 
for the RQ3] and |A04| samples are much more noisy than those 
given in the more recent papers. We note that these two older 
studies are not purely spectroscopic. They both use photometric 
data to estimate effective temperatures and trigonometric paral- 
laxes to determine the surface gravities, therefore they are not 
strictly differential. This is not the case for the more recent pa- 
pers, which use only the spectra to determine the stellar param- 
eters. To make sure that Fig. [3] is not biased by these two older 
stud i es, we re- computed the average abundances but only for 
lT07l IN09L IgToI and HO. The result is shown in Fig.0 where 
it can be seen that the conclusion is not altered by excluding the 



more noise data. The last column of Table [4] gives the average 
abundance ratios plotted in Fig. [4] 

None of the abundance studies scrutinized in this paper dis- 
agree with the M09 and R09 results. In particular, with the ex- 
ception of the R03 results, which we suspect to be in error for 
the volatiles, none of them suggest a significant enhancement of 
volatiles with respect to refractories. Although very good agree- 
ment with our previous results is only found for two of the six 
studies, the other four are still consistent with them. 

The internal errors in the abundance ratios of the six studies 
used in this paper are expected to be in the 0.04 to 0.10 dex 
range , w hich i s similar to the amplitude of the effect detected by 
M09 and R09. The latter, on the other hand, reported abundance 
ratios with a precision of about 0.01 and 0.03 dex, respectively, 
which were achieved thanks to the extremely high quality and 
homogeneity of the data, as well as the high internal consistency 
of the analy sis. Th eref ore, i t is not surprising that the Tc trends 
observed by lM09l and lR09l were more difficult to detect in the 
six studies compiled for this paper. 



4. A signature of terrestrial planet formation? 



As explained by M09 and also IR09L the solar abundance 
anomaly (i.e., the deficiency of refractory elements relative to 
volatiles in the Sun as compared to solar twins) can be inter- 
preted as a signature of the formation of planets, and perhaps 
more specifically terrestrial planets, in the solar system. They 
speculate that the accretion of chemically fractionated material 
into the solar convective envelope during the formation of the 
solar system made its chemical composition pattern depleted in 
refractories as these elements were retained by the planets. The 
fact that the total mass of refractories in the terrestrial planets of 
the solar system today is of the same order of magnitude as the 
mean observed difference between refractories and volatiles in 
the Sun relative to solar twins supports the hypothesis of a ter- 
restrial planet signature, as well as the observed abundance pat- 
tern of meteorites, which mirrors th e solar abundance anomaly 
dAlexander et alJl200UlCieslall2008l) . 

Fig. [3] shows that the abundance ratios are a roughly linear 
increasing function of condensation temperature for Tc <: 900 K 
but they become nearly constant for Tc 900 K. This suggests 
that the process responsible for the observed trend affects only 
elements with Tc <; 900 K0 These high temperatures are only 
found in the inner parts (<2 AU) of protoplanetary disks, regions 
which, at least in the solar system, are inhabited by terrestrial 
planets. Moreover, M09 analyzed the abundance vs. Tc trends 
of about 10 solar analogs with and without detected giant gas 
planets. They found that those stars hosting hot Jupiters all be- 
have like the majority of solar twin stars, i.e., they have a chem- 
ical composition slightly different from solar. More striking is 
the fact that the fraction of stars which closely resemble the Sun 
in their chemical composition is about 70 % if the star does not 
have a close-in giant planet. Within our interpretation of the Tc 
trend as a planet signature, we are led to conclude that the pres- 
ence of hot Jupiters prevents the formation of terrestrial planets. 
One could also speculate that in those systems, the smaller plan- 
ets have already been accreted by the host star, thus removing 
the initially imprinted abundance signature. 



4 Note that our abundance ratios are all relative to Fe, which, within 
our "planet formation signature" scenario, is depleted in the solar pho- 
tosphere. Had we chosen a volatile element as reference, abundance 
ratios of Tc S> 900 K elements would be nearly constant at about zero. 
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4 6 8 10 12 14 

FIP (eV) 

Fig. 5. Condensation temperature (7c) versus first ionization po- 
tential (FIP) for the elements listed in Tables [2] and [3] 



In order for this proposed terrestrial planet signature to be 
imprinted in the chemical composition of the Sun, however, a 
quite short time-scale for the formation of the solar convective 
zone a nd/or an unusually long lived protoplanetary d isk are re- 
quired (lM09llNordlund1l2009l lGustafsson et al.l l2010). Standard 
hydrostatic models for the early evolution of the solar interior 
suggest that it took about 30 Myr for the convective z one to 
reach its present size fe.g.. lD'Antona & Mazzitellll[l9 94) while 
age determinations of observed disks around young stars sug- 
gest that they dissipate in about 10 My r, at most (e.g.. lMamaiekl 
2009: lMever|[2009HFedele et alJl2010h . If this were the case, the 
chemically fractionated material would have been mixed effi- 
ciently in an early massive convective zone and the signature of 
planet formation would have been diluted. Hydrodynamic mod- 
els of early stellar evolution, however, suggest that the Sun was 
never fully convective and that the internal structure of the early 
Sun ( ^ 1 Myr) was not very differ ent from that of the current 
Sun dWuchterl & Ts charnuter 2003). Within this latter scenario, 
the abundance signature can be effectively imprinted, even if 
the protoplanetary disk was short-lived. Another alternative for 
solving this problem involves the inclusion of episodic accre- 
tion in early stellar evolution calculations, which tends to de- 
crease the size of the con vective envelope (Baraffe et al. 2009; 
iBaraffe & Chabrierll2010h . 

The condensation temperatures of elements correlate very 
well with their first ionization potentials (FIP), as shown in Fig. [5] 
for the elements analyzed in this paper. One may therefore ques- 
tion whether the abundance trend with Tq that we attribute to 
planet formation is instead due to an FlP-related effect. The 
abundances of low FIP elements in the (average quiet) upper so- 
lar atmosphere are found to be ab out 4 times gre ater than those 
observed in the photosphere (e.g.. lFeldm an 1992). This so-called 
"FIP effect" is attributed to the acceleration of low FIP ions from 
the chromosphere into higher layers, a process probably driven 
by magnetic fields (e.g., Henoux 1998). The FIP effect is not ex- 
pected to affect the photospheric abundances. Nevertheless, we 
compared the significance of the correlation between [X/Fe] and 
FIP with that corresponding to Tc for the sa mples analyzed in 
this paper as well as those from the M09 and R09 studies. The 
results are given in Fig.|6]and Tables |5]and[6] 

In most cases, the Spearman correlation coefficient (in abso- 
lute value}js larger for the Tq trend than the FIP trend. For the 
|R03| and |A04| samples, these values are similar although both 
correlations (FIP and Tq) are much more significant in thelAOl 




4 6 8 10 12 14 500 1000 1500 
FIP (eV) T c (K) 

Fig. 6. Comparison of the mean abundance ratio [X/Fe] ver- 
sus first ionization potential (FIP) relation (left panels) with the 
[X/Fe] versus condensation temperature (7c) relation (right pan- 
els) for each of the samples studied in thi s paper (top six panels), 
and the Ramirez et al .] d2009l lR09h and iMelendez et ail d2009l 
M09) studies (bottom two panels). 

Table 5. Spearman correlation coefficient r$ and percent prob- 
ability of correlation occurring by chance (P$ ) for the samples 
shown in Fig. [6] 



Sample 


rs 


Ps 


rs 


Ps 




(vs. FIP) 


(vs. r c ) 


R03 


0.125 


54.8 


-0.138 


51.5 


A04 


-0.654 


1.5 


0.648 


1.4 


T07 


-0.090 


55.9 


0.399 


21.1 


N09 


-0.093 


61.9 


0.576 


7.9 


G10 


0.111 


60.7 


0.063 


63.9 


BFO 


-0.164 


50.6 


0.278 


38.2 


R09 


-0.433 


9.7 


0.672 


0.8 


M09 


-0.571 


1.3 


0.736 


0.1 
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Table 6. Spearman correlation coefficient rs , percent probability 
of correlation occurring by chance (P$ ), and significance of the 
[X/Fe] vs. FIP or T c slope for the T c > 900 K data points of the 
samples shown in Fig. [6] 



Sample 


rs 


Ps 


slope ja 


rs 


Ps 


slope/cr 






(vs. FIP) 






(vs. T c ) 




R03 


-0.397 


10.0 


-2.1 


0.380 


13.5 


1.6 


A04 


-0.616 


3.8 


-2.9 


0.626 


2.8 


2.6 


T07 


+0.024 


58.7 


+0.1 


0.386 


25.1 


1.2 


N09 


-0.093 


61.9 


-0.3 


0.576 


7.9 


2.0 


G10 


+0.111 


60.7 


+0.6 


0.063 


63.9 


1.0 


BFO 


-0.154 


55.5 


-0.5 


0.322 


36.3 


0.9 


R09 


-0.222 


46.2 


-1.0 


0.588 


5.1 


2.9 


M09 


-0.370 


17.4 


-1.0 


0.682 


0.6 


2.6 



case. The result for the G10 data suggests a lack of correlation 
with both FIP and Tc- However, this is mainly due to their Sc 
abundances (see below for a quantitative statement), which we 
suspect to have a small but significant error in its zero point 
(Sect. 0. In four of the samples listed in Table [5] the Tc trend 
has a correlation coefficient greater than 0.5. Only two samples 
show similarly high correlation coefficients for the FIP trends. 
In none of the cases investigated here is the correlation with 
FIP much more significant than that with Tc- This situation is 
made even more clear when only the abundances of refractory 
elements are analyzed (Tabled The correlation coefficients for 
the Tc trends suggest a similar or higher significance compared 
to the FIP trends for all samples except G10, unless we exclude 
their Sc abundances that have a suspected small offset, in which 
case rs = 0.054 for FIP and rs = 0.278 for T c ■ F or fou r cases the 
FIP trend is not significant (IT071 HollGlol andfBFO). We also 
note that a linear fit to the abundance vs. Tc trend of refractories 
is, for the majority of cases, more significant than the FIP trend. 
For example, the slope with Tc for the lN09l sample is significant 
at the 2.0 o" level whereas that corresponding to the FIP trend is 
not significant at all. For the G10 results excluding Sc we find 
that the Tc slope divided by its error is about 5 times larger (in 
absolute value) than that corresponding to the FIP trend. Thus, 
we find that, statistically, the [X/Fe] abundance ratios correlate 
bet ter wi th Tc th an FIP. Note that this is particularly true for 
the|MQ9|andR09 results, which are the most reliable due to the 
particularly high abundance accuracy. 

Finally, we should also note that, for all stu dies i n which this 
element is included (the only exception being lA04h . Na, which 
has one of the lowest FIP (5.139 eV) but not a particularly high 
Tc (958 K), fits very well the abundance vs. Tc trend but reveals 
itself as an outlier of the abundance vs. FIP correlation. A similar 
behavior is exhibited by K but, as noted in Sect[3] its abundance 
is not as reliable. 

5. Future directions 

The literature data for solar analogs analyzed in this paper sup- 
port the observational findings of M09 and R09. The uncertain- 
ties in the abu ndan ces of ind ividual stars are significantly larger 
than those by M09 and lR09l and it is therefore not possible to es- 
tablish with confidence which fractio n of s tars have abundances 
extremely close to solar (except for lN09l for which we find a 
value of 10 to 20 %). According to our interpretation, stars with 
no deficiency of refractories relative to volatiles (or overabun- 
dance o f volat iles) relat ive to the Sun should also host terrestrial 
planets. |M09| and |R09| were able to determine a frequency of 
about 15 % for these objects. 



5.1. Metallicity effects 

From studies of extrasolar giant gas planets, it has been well 
established that the fre quency of planet-h osts is roughly propor- 
tional to 10 2[Fe/H] (e.g jMarcv et al.l2005l) . If we assume that the 
relative frequency of planet hosts is independent of whether the 
planets are gas giants or terrestrial planets, the fraction of ter- 
restrial planets at [Fe/H] -0.2 should be about 6% while 
that at [Fe/H] +0.2 it should be about 40 % based on the 
above-mentioned 15 % frequency at [Fe/H] = 0. For a sample 
of about 20 stars, this means that at [Fe/H] -0.2 the number 
of stars which would closely resemble the Sun in their chemical 
composition is 1. At [Fe/H] 0.0 the corresponding number 
is 3, as in IR09L In both cases such a low number makes these 
objects an obvious minority and therefore the solar abundance 
anomaly manifests itself clearly. Following this line of reason- 
ing, at [Fe/H] =s +0.2 the number of stars without Tc trends in 
the abundances should be about 8, a number that is comparable 
to that of those stars which show the Tc trend (12). If the hypoth- 
esis of planet formation signature is correct and the observed fre- 
quency of gas giant planet hosts as a function of [Fe/H] applies 
also to the planets responsible for the abundance signature we 
have uncovered, the Tc trend with (average) abundance ratios 
should progressively disappear at higher [Fe/H]. The signature 
in this case has been washed out by taking the average of a simi- 
lar fraction of stars which show and do not show the Tc trend. If 
sufficient precision is achieved, it may be possible to distinguish 
between metal-rich stars that show and do not show the planet 
signature from a bimodal distribution of Tc slopes for refractory 
elements. 

Interestingly, the behavior of the average abundance trend of 
solar analogs with Tc as a fun ction of [Fe/H] described above 
is in good agreement with the lN09l data (Fig. |7). We perform 
this test on this sample only because it shows the most clear 
and a very significant Tc trend for refractory elements. Fig. [7] 
shows that the solar abundance anomaly is observed both at 
[Fe/H] =a -0.2 and [Fe/H] 0.0. This is expected accord- 
ing to our interpretation because the fraction of stars with abun- 
dances extremely similar to solar (and therefore potential terres- 
trial planet hosts) is very small in both cases. At [Fe/H] +0.2 
the fraction of terrestrial planet hosts is comparable to that of 
stars which do not host those planets and therefore the signa- 
ture is no longer evident. We also expect the scatter to be larger 
in the metal-rich case because of the even mixture of stars that 
have different chemical compositions. Indeed, the weighted 1-cr 
scatter of the average abundance ratios relative to the linear fits 
shown in Fig.|7]increases from about 0.012 dex at [Fe/H] = +0.0 
to 0.024 dex at [Fe/H] = +0.2. We note, however, that the corre- 
sponding scatter at [Fe/H] = -0.2 is larger (0.033 dex) than that 
of the two more metal-rich samples. 

In addition to the result described above. lR09l found that the 
fraction of stars showing the planet signature in the abundances 
increases with [Fe/H] while a hint of a bimodal distribution for 
the Tc slopes of refractory elements was also observed (their 
Fig. 4). Although encouraging, these results rely on the preci- 
sion of abundance ratios relative to solar for stars which are sig- 
nificantly more metal-rich (or more metal-poor in the N09 case) 
than the Sun and therefore small systematic errors and/or chemi- 
cal evolution may affect the abundance determinations. It is cru- 
cial that we repeat the exercise of looking for abundance trends 
with condensation temperature but using large samples (^ 50) 
of metal-rich (or metal-poor) twin stars and analyzing them with 
respect to a well known metal-rich (or metal-poor) reference star. 
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Fig. 7. Upper panel: Average [X/Fe] versus Tc relation for re- 
fractory elem ents in solar analog s with [Fe/H] = -0.2 + 0.1. 
Data are from lNeves et alj ((2009). The dashed line is a linear fit 
to the data while the solid line, whic h corresponds to the aver- 
age relation for solar twins found bv iRamfrez et al.l (2009), has 
been drawn for comparison. Middle panel: as in the upper panel 
for [Fe/H] = +0.0 + 0.1. Lower panel: as in the upper panel for 
[Fe/H] = +0.2 + 0.1. 



5.2. Effective temperature effects 

Our proposed planet signature in the abundances results from 
the accretion and subsequent mixing of chemically fractionated 
material into the stellar convective zone. Thus, the magnitude of 
the effect depends on the size of the convective envelope at the 
time of protoplanetary disk accretion. A massive envelope would 
dilute any chemical peculiarities of the accreted material. It is 
well established that the mass of the convective envelope reduces 
with effective temperature. In F-dwarfs at r e jf 6250 K, the to- 
tal mass of the convective envelope is 0.004 M Q , a figure that is 
significan tly smaller than that cor responding to the Sun (about 
0.02 M ; Pins onneault et al.1 12001). This implies that the chemi- 
cal signature, if due to the accretion from a protoplanetary disk 
of similar mass as the one that gave origin to our solar system, 
should be enhanced by a factor of ^ 5 in mass, everything else 
being equal. This corresponds, on the logarithmic abundance ra- 
tio scale, to a maximum difference in refractory and volatile el- 
ements of 0.30 dex instead of 0.08 dex for the solar case. The 
samples studied in this paper do not exhibit such behavior. Note, 
however, that the analysis of these warmer stars is much more af- 
fected by systematic errors such as those due to departures from 
the local-thermodynamic e quilibrium a ssumption as well as sur- 
face inhomogeneities (e.g.. lAsplundl2005l) . Also, their relatively 
large projected rotational velocities and macroturbulent veloci- 
ties make the identification of line blends more difficult and con- 



sequently the line profile fitting or equivalent width measure- 
ment less accurate. As with the metal-rich example described 
above, it is crucial that we analyze very high quality spectra of 
a large sample 50) of twin F-dwarfs with respect to a well- 
known object of that class in order to search for the abundance 
patterns found for solar twins stars, if indeed they are there. 

5.3. Recent developments 

During the review process of this paper, the two suggestions de- 
scribed above regarding the impact of metallicity and effective 
temp erature on the abundance vs. Tc trends have been tested by 
iGlOl but using a different approach, that of weighting the results 
according to the distance between stars in the stellar parameter 
space. The method that we use is more direct although it remains 
to be proven whether it is more reliable. The size of the samples 
that we propose to analyze (about 50 stars in each case) should 
allow us to make quantitative tests of the planet signature hy- 
pothesis. 

Also during the re view process of this paper, 
iGonzalez Hern andez et al. (2010) published a search for 
our proposed planet signatures in the abundances that they 
derived for a sample of solar analogs. Most of their sample 
stars and data employed are included in lN09i They find good 
agreement between their abundance vs. Tc trends for solar 
twins and the results by M09. They claim, however, that their 
abundance trends for metal-rich solar analogs are incompatible 
with those by R09. In particular, they discuss their results for 
two stars known to host super-Earth like planets (HD 160691 
and HD 1461). They find positive [X/Fe] vs. Tc slopes for 
these two stars, which in our interpretation would suggest the 
absence of rocky planets. In order to understand the apparent 
discrepancy with IR09L in Fig. [8] we show and examine the 
[X/Fe] vs. Tc trend obtained by Gon zalez Hernandez et al.l for 
one of these two super-Earth planet hosts (HD 160691). The 
situation for HD 1461 is very similar, as is that for the majority 
of the other super-solar metallicity stars discussed in that paper. 

The Tc slopes they use to examine our planet formation in- 
terpretation include both refractory and volatile elements. The 
two stars in question are both metal-rich. One must be care- 
ful when interpreting these abundance results beyond the realm 
of solar twins. Chemical evolution is likely affecting the abun- 
dances of the volatile elements C and O much more than the 
majority of refractories at these high metallicities. The other two 
volatile elements in their work are S and Zn, both of wh ich can 
be very sensitive to errors in stellar parameters. In |R09| we ex- 
cluded these elements from the discussion of Tc slopes because 
the trend for refractories is more robust. Note also that, as we 
have pointed out before as being one of the most important char- 
acteristics of the [X/Fe] vs. Tc relation, the overall Tc trend is 
not linear but rather shows a break at Tc - 1000K. 

Fig. [8] shows that the Tc slope for HD 160691 is positive 
if one includes all elements (dashed line)0 As stated above, 
however, the abundances of volatiles are not as relia ble as the 
trend for refractories, which is the reason why |R09| interpre- 
tation is based on Tc slopes for Tc > 900 K elements only. 
Interestingly, a linear fit to the refractory elements only (Tc > 



5 The value that we derive for this slope is +0.057 x lO^dexK 1 , 
which is slightly different from that given by Gonzalez Hernande z et al.1 
(2010) because they mistakenly considered Sr a volatile elem ent, plac- 
ing it a t Tc — 450 K rather than at Tc = 1464 K, as provided bv lLoddersI 
(2003). Note that this is also true for all the T c slopes presented in that 
paper. 
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Fig. 8. Abundance ratios [X/Fe] vs. condensation temperature 
Tc for t he super-Earth planet host HP 160 69 1 . Abundance data 
are from Gonzalez Hernandez et al . (2010). The dashed line is a 
linear fit to all elements; the solid line is a lin ear fit to refracto- 
ries (7c > 900 K) only, as was done by |R09| to investigate the 
planet signature interpretation. Volatiles must be excluded from 
this discussion because their abundances are not as reliable be- 
yond the solar twin regime. 




Fig. 9. Slopes of the [X/Fe] vs. Tc relation for refractory el- 
eme nts only (Tc > 900 K) as d erived by R09 (squares) 
and lGonzalez Hernandez et"aL (1201 Oh (circles). Two super-Earth 
planet hosts are shown with the star symbols. 



900 K) reveals a negative slope of -0.062 x 10~ 3 dexKr' for 
HD 160691 (solid line in Fig. [8} and also a negative slope of 
-0.067 x 10 ~ 3 dexK~' forHD 1461 inst ead of the positive slopes 
obtained by Gonzal ez Hernandez et al.l 

In Fig. [9] we sh ow the distribution of Tr > 900 K slopes for 
both the R09 and lGonzalez Hernandez et"al1 d2010l) samples (we 
d erive d the slopes for the latter using the same procedure used 
in R09 and the abundance data provided in their paper). There is 
good agreement in the distrib ution of Tc > 900 K slopes of these 
two studies even though the iGonzalez Hernandez et al.l sample 
extends to higher metallicities. Moreover, the two super-Earth 
planet hosts, shown with the star symbols in Fig. [9] fall in the 
midst of the distribution of slopes for metal-rich stars. 



The hint of a bimodal distribution of Tc > 900 K slopes 
at hi gh metallicity speculated bv | R09l is not confirmed by 
the IGonzalez Hernandez et al.l (|2010i) data. Instead, beyond 
[Fe/H] =a 0.15, almost all stars show near zero or negative slope. 
Within the planet formation signature in the abundances sce- 
nario, this suggests that all stars with metallicities higher than 
[Fe/H] =i +0.15 are terrestrial planet hosts. In particular, the two 
super-Earth planet hosts HD 160691 and HD 1461 both reveal 
the signature of terrestrial planets in their abundance pattern. We 
thus conclude that, when examined properly, the abundance vs. 
Tc trend s of th ese two super-Earth planet hosts are fully consis- 
tent with |R09| results and interpretation. 



6. Summary 

We have compile d abundance data for six independent spec- 
troscop i c surv eys ( Reddv et a l. 2003; Allende Prieto et all 2004; 
Take dal 120071: " iNeves et al.l l2009t (Gonzalez et all 12010: and 
iBensbv et al.L in preparation) and computed mean abundance ra- 



tios [X/Fe] for solar twins and analogs in each of them. With the 
exception of one study, we find that [X/Fe] correlates with the 
condensation temperature of the elements in a manner similar 
to that found by M09 and R09. We show that it is very likely 
that t he one s tudy (Red dv et al.l 12003) that disagrees with the 
M09 and[R09 results suffers from systematic errors that affect 
mostly the abundances of volatile elements, which are the only 
discrepant ones. Quantitatively, the agreement between the Tc 
trends derived by us in our previous work and those from the 
literature is good for two of the six studies analyzed, while it is 
still consistent for the other four. 

When the abundances from the six studies used in this paper 
are combined without attempting to homogenize the abundance 
results in any way, the agreement with lM09l and lR09l is also 
good. We argue that the heterogeneous nature of this compila- 
tion makes the finding of a good agreement even more robust. 
All six studies are completely independent and employ largely 
different spectra, model atmospheres, spectrum synthesis codes, 
linelists, etc. 

We have investigated whether the abundance trends that we 
find with condensation temperature are due to an underlying 
correlation with the first ionization potential of the elements. 
Statistical tests show that condensation temperature is the driv- 
ing parameter in these correlations. 

The abundance trends with condensation temperature found 
for solar twins by M09 and R09 and shown to be consistent with 
our survey of literature data are speculated to be due to the for- 
mation of terrestrial planets in the Sun and a minority of solar 
twins. We provide suggestions to test specific predictions of this 
hypothesis. 
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